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ABSTRACT 

We report a set of numerical experiments aimed at addressing the applicability of competitive 
accretion to explain the high-mass end of the stellar initial mass function in a sheet geometry with 
shallow gravitational potential, in contrast to most previous simulations which have assumed formation 
in a cluster gravitational potential. Our flat cloud geometry is motivated by models of molecular 
cloud formation due to large-scale flows in the interstellar medium. The experiments consisted of 
SPH simulations of gas accretion onto sink particles formed rapidly from Jeans-unstable dense clumps 
placed randomly in the finite sheet. These simplifications allow us to study accretion with a minimum 
of free parameters, and to develop better statistics on the resulting mass spectra. We considered both 
clumps of equal mass and gaussian distributions of masses, and either uniform or spatially-varying 
gas densities. In all cases, the sink mass function develops a power law tail at high masses, with 
dN/dlogM oc M~ T . The accretion rates of individual sinks follow M oc M 2 at high masses; this 
results in a continual flattening of the slope of the mass function towards an asymptotic form r ~ 1 
(where the Salpctcr slope is T = 1.35). The asymptotic limit is most rapidly reached when starting 
from a relatively broad distribution of initial sink masses. In general the resulting upper mass slope is 
correlated with the maximum sink mass; higher sink masses are found in simulations with flatter upper 
mass slopes. Although these simulations are of a highly idealized situation, the results suggest that 
competitive accretion may be relevant in a wider variety of environments than previously considered, 
and in particular that the upper mass distribution may generally evolve towards a limiting value of 

r~i. 

Subject headings: stars: formation — stars: luminosity function, mass function — ISM: clouds 



1. INTRODUCTION 

The stellar initial mass function (IMF) among other 
things determines the fraction of stellar populations in 
massive stars; this in turn affects the production of heavy 
elements, the stellar feedba ck of energy into the ISM, and 
the evolution of galaxies. ISalpeterl ( 1955|) first pointed 
out the power law distribution in the "original mass 
function" ; subsequent observational work has established 
the general form of the IMF, which at high masses 
is still comparable to the "Salpeter slope" T, where 
dN/dlogM = M~ r , T = 1.35. The most widely used 
functional form is a power-law distribution or a combi- 
nation of power-law distribution at different mass ranges. 
Other widely used forms of the IMF include log-normal 
distributions and combi nation of power-law and log- 
normal di stribution (e.g . IChabrieil l2003t IBastian et al.l 
2010). As IBonnell etaLl (1200717 noted, the essential fea- 
tures of the IMF include a peak at a mass of a few tenths 
of Mq and a declining power-law tail toward higher 
masses. 

While the origin of the IMF remains a matter of ex- 
tensive debate, two general i deas have com e to promi- 
nence in recent years (e.g., iClarkd I2009D . The first 
supposes that the mass spectrum of dense structures 
within star-forming clouds, suggested to be the result 
of supersonic turbulence, more or less directly maps into 
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the s t ellar mass d i stribu tion (e.g., iPadoan fc Nordlundl 
I2002t iKlein et al.l I2007D . In these models the IMF 
results fr om local mass reser v oirs that are relatively 



isolated (jPadoan et al.l I2007t IHennebelle fc Chabrierl 

la 



2008 ), possibly affected by gravity (jKlessen et al 
iKlessen &: Burkertll200lT ). The second type of model in- 
vokes two processes to produce the IMF; the low-mass 
end is determined by turbulence and thermal physics, 
qualitatively similar to the first picture, but the high- 
mass "tail" is a result of continuing accretion from a ma ss 
reservoir (e.g., IZinneckerl [1981 IBonnell et all I2001allb1 . 
I2007h . Thus the accumulation of material by the most 
massive stars is the result of non-isolated accretion, from 
size scales greater than the local Jeans length. The pro- 
cess resulting in producing the high-mass end of the IMF 
in this approach is usually called "competitive accretion" 
(CA). 

As summ a rized by Cla rk et "ail (1200 !); and 
IBonnell et all (|2007f ). the high -mass power-law tail 
in CA simulations typically arises from formation in 
a stellar cluster; the potential well results in high 
gas densities near the center, helping to f eed materia l 
into the most massive objects (see also iBatd [2009). 
Bonnell et al. (2001b) found that the slope of the 
mass function depended upon whether the gravitational 
potential was dominated by gas - in which case they 
found an asymptotic limit of T = 0.5, due to tidal lobe 
limitation of mass accretion; or by stars, in which case 
the asymptotic limit was r = 1, where Bondi-Hoyle 
accretion dominates. The latter is consistent with the 
analysis of Zinnecker (1982), who showed that T = 1 
results asympotically from accretion rates which scale 
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as M oc M 2 . 

These investigations suggest that CA can account for 
the high-mass end of the IMF in clusters. However, while 
most stars form in clusters, a non-negligible number do 
not, at least in the solar neighborhood. In addition, the 
properties of clusters vary widely, with most being rela- 
tively small (Lada & Lada 2003); this raises the question 
as to whether the IMF might be affected by the mass 
of the cluster. Moreover, the initial states and evolu- 
tion of protocluster clouds and clusters are uncertain; 
current assumptions range from relatively slow evolution 
in a roughly virialized condition (e.g., Tan, Krumholz & 
McKee 2006) to the the opposite assumption of rapid 
gravitational collapse (e.g., Tobin et al. 2009; Proszkow 
et al. 2009). We are therefore motivated to investigate 
a schematic model of competitive accretion which does 
not employ the assumption of formation in an initially 
clustered environment. In addition, we wish to adopt 
a simple initial physical model with as few parameters 
as possible to isolate the most important properties for 
producing the high-mass IMF. 

In this paper we report a set of numerical simulations 
in a simplified model to address some general aspects of 
competitive accretion. Our results suggest that values 
of T close to the Salpeter slope can result in a wider 
variety of environments than previously discussed; they 
also suggest that the value of V may be correlated with 
the maximum mass achieved through CA. These findings 
suggest additional new approaches for numerical simula- 
tions of the production of stellar IMFs. 

2. MODEL AND METHODS 

Our initial setup is motivated by our models of 
molecular cloud formation as a result of large-scale 
flows in the interstellar mediu m dHeitsch et al.l 120061: 
Vazquez-Semadeni et al.l [20061 : iHeitsch fc Hartmannl 
20081 iHeitsch et al.ll2008allbT r In these models the dense 
material formed in post-shock gas is geometrically thin 
rather than spherical, due to post-shock compression by 
large-scale flows (e.g., Hartmann, Ballesteros-Paredes 
& Bergin 2001). As there is no particular mechanism 
which would enforce virialization, the cloud as a whole 
collapses laterally under gravity; eventually, much if 
not most of the supersonic motion in the cloud is due 
to acceleration by the cloud's self-gravity, rather than 
the initial turbulent velocities injected during cloud 
formation (e.g., Heitsch et al. 2008; Heitsch & Hartmann 
2008). The most important role of this mostly gravita- 
tionally driven turbulence in the post-shock gas is to 
provide density enhancements which can gravitationally 
collapse faster than the cloud as a whole (Heitsch, 
Hartmann, & Burkert 2008). 

We adopt an extremely simplified version of this cloud 
formation model; specifically, we use an initially circular 
isothermal sheet with many thermal Jeans masses ini- 
tially in hydrostatic equilibrium in the short dimension. 
We then introduce local Jeans-unstable mass concentra- 
tions in a spatially-random pattern within a given ra- 
dius which rapidly form sink particles (protostars). For 
simplicity we do not introduce initial velocity perturba- 
tions; instead, we allow the cloud and sinks to evolve 
under their own gravity. The random placement of the 
sinks (along with any density fluctuations imposed in the 
gas) quickly results in complex "turbulent" gas velocities 



which are gravitationally-generated. This setup allows 
us to avoid the issue of fragmentation for the present 
and concentrate on the development of CA in an ini- 
tially non-clustered environment with a minimum of free 
parameters. 

We use Gadget-2 (jSnringel et al]l2"0?Hl : ISpringelll2T)05h 
to simulate the gas dynam ics and the formation of "pro- 
tostellar" sink particles. Uappsen et al.l (|2005l ) imple- 
mented the sink particle formulation into the form of 
Gadget-2 we use. Collapsing structures above a density 
threshold (n = 10 cm -3 in our case) are replaced by sink 
particles, which interact with gas and other sink particles 
through only gravity. 

For simplicity, we assume an isothermal equation of 
state at 10 K for the gas particles, with a molecular 
weight of [a = 2.36. We use a code unit system in which 
the unit length is 1 pc, the unit time is 1 Myr and the 
unit mass is 0.058 M . In these units the radius of the 
sheet is then 2 pc and the total mass of the sheet is 
820 M Q . The surface density of the unperturbed sheet is 
1.37 x 10 _2 gcm -2 (Ay — 3.8 perpendicular to the sheet). 
The (initial) number of gas particles in each simulation is 
Ntot = 1.6 x 10 6 . For convenience we report results scaled 
to the above physical units, but note that the simulations 
can be rescaled given the assumed isothermal equation 
of state. Specifically, if the unit length is scaled to d pc, 
the unit of time becomes d Myr and the unit of mass 
becomes 0.058 dM Q . 

The initial vertical structure of the sheet follows 

p(z) = p scch 2 (z/H), (1) 

with po — 3.7 x 10~ 20 g cm" 3 and scale height H = 0.06 
pc. However, the equilibrium density distribution of an 
isothermal infinite sheet will follow the same form, with 
a scale height of H = c^ttGS)" 1 = 0.04 pc. 

In x and y directions, the gas particles are randomly 
placed in a uniform sheet, with a radius of 2 pc (except 
for the non-uniform sheet case, see Section 3.2). This 
leads to density fluctuations due to the random position- 
ing of the particles. To plot the surface density and ve- 
locity fields, we interpolated the densities and velocities 
of the SPH particles onto a rectangular grid. Each cell 
has an area of (0.015) 2 in code unit or (0.015pc) 2 . 

We start each simulation with 100 Jeans unstable 
clumps. The rapid collapse of these clumps leads to dy- 
namic creation of sink particles before 0.1 Myr. We were 
unable to put sinks in at the start, probably because of 
problems with the boundary conditions around the sinks; 
when the sinks are dynamically created within the sim- 
ulation, the boundary conditions are properly calculated 
to account for the disco ntinuities in density and gas pres- 
sure around the sinks (jBate et al.l fl995: Jap psen et al.l 
120051) . 

Because the sheet itself is also highly Jeans unsta- 
ble, it also collapses under gravity, on a timescale t c ~ 
RinGY.)- 1 / 2 w lAMyr (jBurkert fc Hartmannl l2CJ0l 
hereafter BH04). Due to gravitational focusing, a ring of 
material piles up quickly along the edge of the cloud. The 
edge can then become gravitationally unstab le and frag- 
ment fBH04: IVazquez-Semadeni et all 120071 : Figure [TJ. 
With our isothermal equation of state, we find relatively 
uncontrolled (numerically) fragmentation in this ring; we 
therefore turn off the creation of sinks after the initial 100 
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clumps collapse, allowing us to focus entirely on compet- 
itive accretion within the main body of the cloud. Our 
restriction on the initial placement of clumps to a radius 
of 1 pc avoids accretion from the ring. 

Gas particles that come within a certain radius of a 
sink (0.003 pc in our setup) are tested for accretion in- 
dividually. If a gas particle is bound to a sink, the gas 
particle is accreted by the sink. Gas particles which come 
within 0.0003 pc of the sink are always accreted. We ran 
each simulation for 1.2 Myr, or approximately 0.8 t c , with 
an output file written every 0.1 Myr. 

Within this general setup we considered several cases. 
In the first set of simulations, we assumed a uniform 
surface density for the cloud and that each clump had 
the same mass, 0.82M©. In a second set, we assumed 
the same equal initial clump masses, but a varying den- 
sity distribution in the gas. The final sets of simulations 
assumed constant surface density gas but log-normal ini- 
tial mass distributions for the clumps, keeping the total 
mass of the clumps to be 10% of the cloud mass. To 
improve statistics, we ran six realizations of each of the 
simulations described above, differing only in the random 
positions of the clumps. 

3. RESULTS 
3.1. Equal Mass Clumps in a Uniform Sheet 

Figure Q] shows one of the realizations of the simplest 
case, equal mass clumps in a uniform sheet. The left 
panel shows the view from the top, and the right panel 
shows the side view. Figure [2] shows a close-up view of 
the central 1.2 x 1.2 pc. The circles mark the location of 
the sink particles, and the area of the circles correspond 
to the mass of the sinks. 

Early on (before 0.2 Myr), most sinks evolve indepen- 
dently of each other, accreting mass from the original 
clump and the environment. However, as the entire cloud 
collapses, after 0.2 Myr, the sink particles start to affect 
each other, forming small groups, in a manner reminis- 
cent of the simulations of Bonnell, Bate, & Vine (2003) 
(see also Maschberger et at 2010). By 0.5 Myr, the gas 
between the sink particles starts to form a filamentary 
structure that resembles the "cosmic web" in cosmologi- 
cal simulations. At this stage, part of the gas is accreted 
first onto the filament, and then from the filament to the 
sinks. The regions between the web become depleted of 
gas. As time goes on, the small groups collapse, creating 
larger groups while the sink particles accrete gas from 
the environment. The more massive sinks in a group can 
accrete mass faster, thus broadening the mass distribu- 
tion 

Figure [3] shows the growth of each sink particle as a 
function of time. Initially, all the clumps have the same 
mass, but the final sink masses span over 1.5 dex in mass. 
Note that the initial clump mass is not equal to the sink 
mass when the sinks are created because it takes about 
0.2 to 0.4 Myr for the all the clump gas to fall in. 

Figure 0] shows the mass accretion rate of each sink at 
intervals of 0.2 Myr, including the sink particles from all 
six runs. The accretion rates of the more massive sinks 
exhibit a roughly dM/dt oc M(sink) 2 behavior. As the 
system evolves, the accretion rates decrease due mostly 
to the removal of gas into sinks, and the lower-mass sinks 
lose the competition for material to the high-mass sinks. 



As shown in Figure [SJ in an initially non-clustered en- 
vironment, the accretion rate shows no clear dependence 
on the position of the sink within the sheet. This is 
unsurprising given the uniform nature of the sheet, al- 
though the global motions of the sheet do depend upon 
radius. This is in contrast to formation in an initially 
clustered environment, as described by Bonnell et al. 
(2001), where the accretion rate depends on the position 
of the sink in the cluster through the tidal lobe radius. 
While the center of the cluster is the preferred location 
to form the most massive star, the most massive stars 
in our simulations do not necessarily form in the center 
(though eventually everything collapses to the center). 

The combined mass distribution of the six runs is 
shown in Figure [6j The thin black line represents the 
initial mass of the clumps. The thick lines show the 
mass distribution at 0.4, 0.6, 0.8, 1.0 and 1.2 Myrs after 
the beginning of the simulation. The distribution starts 
with a delta function, evolves into a Gaussian-like distri- 
bution and then develops a high-mass power-law toward 
the end of the simulation. The solid black line show a 
fit to the distribution from 10 4 to 10 1 " 15 M Q when we 
terminate the simulation, or at t — 1.2 Myr, with a slope 
of — 2.07±0.15. The derived slope does depend modestly 
on the range of masses which are fitted. The slope and 
the fitting range in mass is tabulated in Table [T] 

3.2. Equal Mass Clumps in a Non-uniform Sheet 

The setup is mostly the same as the previous case, 
but with background density fluctuations. To construct 
a varying surface density, we used the linear superposi- 
tion of sine waves in both the x and y directions whose 
magnitude is proportional to the wavelength: 

d(x,y) = 2J k^ 1 sm(k x x + 4> x (k x ))sin(k y y + (f> y (ky)), 

where d(x, y) is the surface density at location x, y; k x 
and k y are the wavenumbers in x and y directions; ip x 
and 4> y are the randomly chosen phases. The /c _1 fac- 
tor is used simply to ensure that the fluctuations are 
mostly on large scales while still having noticeable ef- 
fects on smaller scales. On the smallest scales, the den- 
sity fluctuations are dominated by random positioning 
of the particles. The largest wavelength allowed is the 
diameter of the sheet; the smallest wavelength allowed is 
1/20 of the diameter. The fluctuating part of the surface 
density is then added to a constant surface density part 
so that the minimum density is 30% of the maximum 
density. The phases of the surface density are randomly 
chosen for each of the six simulations. Figure [7] shows a 
close-up view of the central 1.2 x 1.2 pc of one of the runs 
of this case. The fluctuations in the background density 
are not very prominent in the figure partly because the 
surface density is plotted on a log scale, and the clumps 
are dominating the density fluctuations. 

In this set of simulations, the accretion rate is again 
proportional to M (sink) 2 for the more massive sinks (Fig- 
ured]). The sink mass distribution grows in a similar way 
as in the previous case, but the distribution spreads to 
higher masses slightly faster. At t — 1.2 Myr, the linear 
fit to the distribution gives a slope of — 1.97± 0.15, with 
a fitting range of 10 4 to 10 115 M . Thus, including 
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these density fluctuations in the simulation makes little 
difference to the final result. 

3.3. Clumps with an Initial Mass Distribution 

The previous results suggested that a wider initial dis- 
tribution of masses should grow the power-law tail faster. 
We therefore constructed three sets of simulations with 
initial mass distributions 



A(log M) oc exp 



(logM -logA/ c 
2a 2 



(2) 



where log(M c /M Q ) = -0.1 and a = 0.05, 0.1 and 0.2 
dex. Figure IH1 shows the sink mass distributions for these 
three cases. The thin black lines represent the initial 
clump mass distributions, and the thick lines are the 
mass distribution at t=0.4, 0.6, 0.8, 1.0 and 1.2 Myr. The 
wider distribution of initial clump masses yields faster 
growth of the high mass power law, as expected. Lin- 
ear fits to the final mass distribution at t=1.2 Myr yield 
slope of -1.67±0.15, -1.42±0.14 and -1.03±0. 16.' Again, 
the parameters for the fitting are tabulated in Table [TJ 

The slope of the mass distribution depends on the 
spread of the initial clump masses. The final slope can 
be flatter than the Salpeter value of -1.35. In fact, if 
the mass accretion rate grows strictly as M oc M 2 , all 
the slopes would approach -1 if the si nks have enoug h 
time and enough gas to accrete (e.g. , IZinneckeil 119821 ) . 
Our numerical results are consistent with an asymptotic 
slope of r = 1.0, although the statistical errors are large 
enough to prevent an absolutely secure conclusion, even 
with simulations totalling 600 objects. This emphasizes 
the long-standing problem of achieving sufficient num- 
bers of objects, either theoretically or observationally, to 
make firm statistical conclusions about IMF slopes. 

4. DISCUSSION 
4.1. Accretion and clustered environments 
Figures 2] and [5] show the main result of this paper: 
a strong tendency for M oc M 2 to develop at the high- 
mass end of the sink mass distribution, in initially non- 
clustered, flat, collapsing cloud environments. This re- 
sults in a general tendency for the high-mass power low 
to approach r = 1 asymptotically, depending upon how 
much mass the sinks can accrete beyond their initial val- 
ues, as shown in Figures [6] and HI To put this in context, 
we constructed a simple analytic model where an initial 
Gaussian distribution of masses is modified by accretion 
with M oc M 2 = aM 2 , where a is a constant. For an 
initial mass Mo, the mass grows as a function of time 



M(t) = 



Mo 



1 - aM Q t 



(3) 



(| Zinneckerl 1 1 9 8 21 ) . The resulting mass grows as M(t) — > 
oo as t — > too — (aM ) _1 . Figure [TU] shows how the 
mass distribution grow with time, plotted in increments 
of 0.16f oo = (aMo)" 1 . This does a suprisingly good job 
of reproducing the numerical simulation results, if accre- 
tion is stopped at differing times. Even though the sim- 
ulation accretion rates do not scale exactly as M 2 , with 
the lower-mass sinks accreting more slowly, this makes 
little difference on the resulting mass distribution. This 
comparison emphasizes that the "competitive" effect in 



CA is not only starving the low-mass systems at the ex- 
pense of the high-mass objects; in terms of producing 
the high-mass power-law, it is the result of differential 
accretion, enhancing the rates at which the higher-mass 
sinks accrete. 

While our starting conditions do not assume an initial 
clustered structure or a deep central gravitational poten- 
tial, our assumed cloud symmetry and lack of turbulence 
or rotation results in forming a cluster of sinks at the 
center. However, the high-mass tail of the mass func- 
tion is strongly developing well before the final central 
cluster is formed. Indeed, we observe M oc M 2 at the 
earliest stages in our simulations, where the clustering is 
minimal (we also see this in a simulation with sinks in 
a uniform sphere - unsurprisingly). It does appear that 
some local grouping is necessary to achieve enough dif- 
ferential accretion to develop a clearly asymmetric mass 
function, based on simulations (not presented here) that 
show when the sinks are initially placed further apart, 
the groups take longer to form and the high-mass tail of 
the IMF evolves more slowly. 

In our simulations, the local groupings happen 
relatively qu i ckly c ompared to the simulation of 
iBonnell et al.l ()2001bD . This is probably because the re- 
laxation time in a sheet be faster than in a sphere of th e 
same central density and total mass (e.g. . IRvbickil[l97lh . 

4.2. Applicability of Bondi-Hoyle accretion 

From their simulations of formation in a cluster po- 
tential, Bonnell et al. (2001a,b) argued that there are 
two regimes of accretion. The first phase was where the 
gravitational potential of the cluster gas dominated, and 
accretion was tidally limited, leading to a T ~ —0.5. 
This occurs when both the protostars and the gas both 
fall in toward the cluster center (see, e.g., discussion in 
Clark et al. 2009, §2). During the second phase, the 
stars dominate the potential, become virialized, and then 
Bondi-Hoyle accretion leads to an upper mass distribu- 
tion r i. 

In contrast, we find r — > 1 even during global collapse, 
for a situation where the infall velocities tend to be larger 
at large radii and the average density is roughly con- 
stant with position (see also Burkert & Hartmann 2004) . 
This occurs as the groups begin to dominate the local 
gravitational potential and generate significant relative 
velocities of the sinks and the infalling gas. This may 
provide local environments equivalent to the global sec- 
ond accretion regime of Bonnell et al. (2001b). The tidal 
limiting phase is much less important in our simulation 
because of the shallower gravitational potential gradient 
of the sheet, so that the characteristic Bondi-Hoyle ra- 
dius of accretion (see below) is always smaller than the 
tidal radius. 

In the simple, isolated version of Bondi-Hoyle accretion 
in three dimensions, 



M ex P R 2 ac 



(4) 



where p is the gas density and v is the (assumed super- 
sonic) relative velocity of the gas and sink, both averaged 
at the accretion radius 



R acc oc GM/v 2 



(5) 
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This results in the usual scaling 

M oc M 2 pv~ 3 . (6) 

Initially, we thought that in our adopted flat geometry 
the accretion rates might scale as 



M oc 



(7) 



where £ is the gas surface density of the sheet; this would 
imply 

M oc MSir 1 . (8) 

In fact, the accretion of the sink particles is more like a 
3D than a 2D flow. This is because the accretion radius 
is effectively embedded in the sheet. In the small groups, 
the velocity dispersion amongst the sinks is about I- 2 
kms -1 . The accretion radius is then 



Rac 



0.08 



/ M 



V 10M G 



I km s 



pc. 



(9) 



From the above equation, we conclude that for sink 
masses up to IOMq, the accretion radius is in gen- 
eral smaller than the scale height of the sheet. Thus 
the mass flow is (non -spherical) Bondi-Hoyle accretion 
(|Bondi fc HovldU94l . 

It is worth noting that our sheets are undoubtedly 
much thinner than realistic molecular clouds. Thus, our 
results suggest that formation of clouds by large scale 
flows, which tend to produce flattened clouds (see §4.3), 
does not alter the basic applicability of Bondi-Hoyle ac- 
cretion for the upper mass IMF (though conceivably the 
results might be different in filament geometry). 

It is difficult to apply the standard formula ((6]) to 
our numerical results because the background medium 
rapidly becomes strongly perturbed. The gas motions 
are not uncorrelated with the sink velocity, as assumed 
in the development leading to equation (|6]), but instead 
tend to be focused toward mass concentrations. The lo- 
cal gas density distribution is also highly perturbed, with 
strong, gravitationally-accelerated flows into and along 
filaments. Bonnell et al. (2001b) attempted to deal with 
these difficulties through the following argument. Con- 
sider a point mass at radius R in some environment, with 
infall velocities 

Vrel OC R~ V (10) 

and gas densities 



p oc R 



(11) 



With these assumptions Bonnell et al. found 

M oc g(t)M 2 R 3 ^ , 



(12) 

where g(t) is a function which allows for the assumed 
homologous evolution of the cluster. This analysis results 
in r — > 1 for sinks whose masses are initially uncorrelated 
with position; Bonnell et al. (2001b) suggested that the 
slope might be steeper if the higher-mass objects reside 
preferentially in the cluster center. 

To see whether the densities and velocities correlate 
with sink mass, we evaluate these quantities at two 
radii: first, at a radius of 2GM/c 2 s , the maximum ac- 
cretion radius in the Bondi accretion formulation in the 
case where the relative velocity between the sink and 
the gas is subsonic; the other at a radius of 0.024 pc, 



which is the distance sound waves can travel in 0.1 Myr 
(the time between snapshots). Figure [TT] shows scat- 
ter plots of sink masses vs. velocities relative to the 
gas, gas density, surface density and surface density di- 
vided by v, with all properties evaluated at R=2GM/c 2 
at t=0.6Myr. Figure [12] shows the same plots, with 
gas properties evaluated at 0.024pc away from the sink. 
The results show that the densities and velocities of the 
gas are not strongly correlated with the individual sink 
masses. Therefore, the accretion rate scales as M oc M 2 . 
This may be a result of having a group of accreting sinks 
experiencing the same environment, as in the discussion 
leading to equation [12] whatever sets the local density 
and flow velocity, the capture cross-section will still scale 
as M 2 . 

This suggests that the important factor is not the 
form of the initial density and velocity distribution but 
whether the global features are uncorrelated with the in- 
dividual sink masses, as in equation (|12[) . In this view as 
long as a group of objects of differing mass "see" the same 
conditions- gas densities and velocities - their differential 
accretion rates will scale as M 2 (the proportionality due 
to the gravitational cross-section). This only holds for 
the most massive objects in each group; the low-mass 
sinks are starved of material to accrete. More gener- 
ally, the absolute value of the mass accretion rate may 
vary from group to group; but as long as each group 
can set up an M oc M 2 relative accretion rate with dif- 
fering constants of proportionality, one may argue that 
the summed population will still asymptotically evolve 
towards T = 1. 

4.3. Turbulence 

Most simulations of star-forming clouds invoke an 
imposed turbulent velocity field, in view of the su- 
personic spectral line widths observed in molecu- 
lar tracers. Our simplified approach, in which 
we do not impose initial velocity fluctuations but 
initial density perturbations, is motivated by re- 
cent simulations which form t urbulent star -f ormin g 
clouds from large-s c ale fl o ws iHeitsch et al.1 12 006) ; 
Heitsch fe Hartmannl (12 008); Hc itsch et al l (|2008bD and 
Vazcmez-Semadeni et al.l (]2006l 120071 ) . These simula- 



tions found that while hydrodynamically-generated tur- 
bulence in the post-shock gas dominates the cloud struc- 
ture and motions in early phases, gravitational ac- 
celeration dominate s the motions at late stages (e.g., 
IHeitsch et al.ll20~08a|) . Similar behavior is seen in models 
in which the turbule nce is not continu ally driven but al- 
lowed to decay (e.g. JBate et al.ll2003D . Thus, the initial 
turbulence provides density fluctuations or "seeds" which 
then generate supersonic motions as a result of gravita- 
tional forces in clouds with many thermal Jeans masses. 
Our models take this view to a simple extreme, where we 
let gravity do all of the (supersonic) acceleration of the 
gas given initial density fluctuations (our clumps). 

The assumption that the largest "turbulent" motions 
are mostly gravitationally-driven is an essential part of 
the competitive accretion picture. Krumholtz et al. 
(2005) argued that the supersonic velocity dispersions of 
molecular clouds are too large for Bondi (and thus com- 
petitive) accretion to be effective; however, this assumes 
that the "turbulent" motions persist and are spatially un- 
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correlated with the accreting masses. In contrast, even 
though large (and roughly virial) velocities develop in 
our simulations, competitive accretion still operates be- 
cause the motions are largely the result of gravitational 
infall to groups, plus global, spatially-correlated collapse 
of both the sheet gas and the sinks. These considerations 
emphasize the importance of understanding the nature of 
"turbulence" in star-forming clouds. 

4.4. Mass functions 

Recently, there have been suggestions that the stel- 
lar IMF is not universal; in particular, that the most 
massive star in a region depe nds upon its richness 
( ([Kroupa fe Weidner|[2003ll2005[ ); also Weidner, Kroupa, 
& Bonncll 2010 and references therein). The models pre- 
sented here also result in a non-universal upper-mass 
IMF, with a suggestion that r ~ 1 is an asymptotic 
limit which is approached most closely when the matter 
accreted is much larger than the initial "seed" mass; and 
thus, to some extent, the slope may correlate with the 
most massive object formed. This is difficult to ascertain 
observationally, in part because of the tradeoff between 
upper mass slope and truncation mass (e.g. , Maschberger 
fc Kro upa 2009). Using the simulations of B onnell et al.1 
(2003) and Bonnell, Clark, & Bate (2008), Maschberger 
et al. (2010) found global values of T slightly greater than 
unity, and V ~ 0.8 in the richest subclusters. This may 
be consistent with our findings of a correlation between 
slope and upper mass. 

It may be worth noting two other situations in which 
r ~ 1 mass functions are found: dark-matter halo sim- 
ulations (below the upper- mass cutoff; e.g., Jenkins et 
al. 2001); and star cluster mass distributions (Elmegreen 
& Efremov 1997; McKee & Williams 1997; Zhang & Fall 
1997; Chandar 2009), although some estimates yield flat- 
ter power- law slopes (e.g., Maschberger & Kroupa 2009). 
Gravitational accretion thus could potentially provide a 
unified explanation of the similarities in these mass func- 
tions. 

5. CONCLUSIONS 

This paper presents numerical experiments using SPH 
simulations to address the general applicability of com- 
petitive accretion in initially non-clustered environments. 
A flat geometry is used to construct a shallow gravi- 
tational potential as opposed to the spherical clustered 
potential used in previous simulations by Bonnell et al. 
(2001a, b). The simplified setup consists of only the most 
important elements in forming the high-mass IMF: differ- 
ential gas accretion onto protostars under gravity. With 
this setup, we were able to produce the high-mass end 
of the IMF with slopes comparable to the Salpeter slope 
r = 1.35. The simple setup also allows us to under- 
stand the mass growth of sinks in details without worry- 
ing about fragmentation and thermal physics, and also 
permits us to generate reasonably statistically-significant 
results for upper mass function slopes. 

The mass growth rate of the sinks follows M cx M 2 
for all high mass sinks, while low mass sinks sometimes 
accrete at lower rates. The high-mass end of the IMF de- 
velops a power-law tail and flattens, with an asymptotic 
slope of r = 1. Variations in initial clumps masses and 
surface density help the power-law tail to flatten faster. 



In our simulations, most systems do not reach the asymp- 
totic slope due to gas depletion. In real molecular clouds, 
stellar feedback as well as gas depletion can terminate the 
accretion and determine the final high-mass IMF slope. 

The present set of simulations are obviously quite ide- 
alized. Our purpose was to elucidate the basic physics of 
CA in as easily-visualized and interpretable a situation 
as possible. The next steps, which are currently under 
way, are to start with more complex density distribu- 
tions and allow sink formation and consequent evolution 
in more complex geometries, and include velocity fields as 
necessary. While we suspect that the physics of compet- 
itive accretion will remain the most important factor in 
creating the high-mass region of the IMF, as previously 
argued by Bonnell et al. (2001a, b, 2003), and Clark et 
al. (2009), further study is needed. 

We wish to thank the anonymous referee for very help- 
ful comments which improved the paper substantially. 
We thank Jeremy Hallum for his efforts in maintaining 
the cluster on which these simulations were computed. 
This work was supported in part by NSF grant AST- 
0807305 and by the University of Michigan. 
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2 
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gaussian(<7 = 0.1) 
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5 
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gaussian(cr = 0.2) 
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-1.03 ± 0.16 
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Figure 2. The same simulation as in Figure[T] but only the inner part is shown. Each box is 2.4 pc by 2.4 pc. The arrows indicate the 
velocity vectors of the gas, with lkrns" 1 marked on the upper right corner of each panel. The colors correspond to the log of column 
density in g cm -2 . 
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Figure 3. The evolution of sink mass as a function of time in one of the runs of the equal mass case. 
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Figure 5. The accretion rate vs. distance of sink from the center of the sheet at t=0.2, 0.4, 0.6, 0.8, 1.0 and 1.2 My in the runs with 
equal mass sinks and the uniform surface density cloud. All 600 sinks from the six runs are included. 
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Figure 6. Left: The distribution of mass at t=0.4, 0.6, 0.8, 1.0 and 1.2 Myr for the constant background density, equal clump mass case. 
The weighted linear fit has a slope of -2.08. Right: The distribution of mass at t = 0.4, 0.6, 0.8, 1.0 and 1.2 Myr for the varying background 
density, equal clump mass case. The weighted linear fit has a slope of — T = -1.95. The dotted lines represent the Salpeter slope — T = 
-1.35 (see Table 1). 
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Figure 7. The top view of the inner part (each box is 2.4 by 2.4 pc) of the case with equal clump mass and background density fluctuation. 
The arrows indicate the velocity vectors of the gas, with lkms - 1 marked on the upper right corner of each panel. The colors correspond 
to the log of column density in g cm -2 . 
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Figure 8. The accretion rate vs. mass of sinks at 0.2, 0.4, 0.6, 0.8 and 1.0 Myr in one of the runs of the fluctuating background case. The 
accretion rates for the higher- mass sinks follow M <x M 2 ; the surface density fluctuations have little effect (see Figure 4 for comparison). 
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Figure 9. Resulting mass distributions for initial clump masses with Gaussian distributions (Top left: a = 0.05; top right: a = 0.1; 
bottom left: a = 0.2) with constant background density. The linear fit slopes are —Y = -1.66, -1.42 and -1.03, respectively. The thin solid 
lines represent the initial clump mass distribution. The dotted lines represent the Salpeter slope — Y = -1.35. 
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Figure 10. The mass distribution at different time for M = aM 2 . The top figure has a narrower initial mass range (a = 0.05 dex) than 
the bottom figure(<r = O.f dex). The thin solid line represents the initial distribution of masses, and the colored lines represent the mass 
distribution in increments of 0.16too = (aMo) . The dotted line represent the Salpeter slope —F = -1 .35. 
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Figure 11. Mass of sinks vs. the gas properties around the sinks for all the cases in the equal initial mass, uniform density case at t = 
0.6Myr. The gas properties are evaluated at Race = 2GM/c^ from each sink. Top left: mass vs. gas velocity relative to the sink, top right: 
mass vs. gas density, middle left: mass vs. surface density, middle right: mass vs. surface density/velocity, bottom left: mass vs. p/v 3 . 
There is no obvious correlation between the gas properties and the sink mass. 
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Figure 12. Mass of sinks vs. the gas properties around the sinks for all the cases in the equal initial mass, uniform density case at t = 
0.6Myr. The gas properties arc evaluated at R = c s X O.lMyr = 0.024 pc from each sink. Top left: mass vs. gas velocity relative to the 
sink, top right: mass vs. gas density, middle left: mass vs. surface density, middle right: mass vs. surface density/velocity, bottom left: 
mass vs. p/i> 3 . There is no obvious correlation between the gas properties and the sink mass. 



